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Abstract 

We show that the viscous Burgers equation ut + uu x = u xx consid- 
ered for complex valued functions u develops finite-time singularities 
from compactly supported smooth data. By means of the Cole-Hopf 
transformation, the singularities of u are related to zeros of complex- 
valued solutions v of the heat equation vt = v xx . We prove that such 
zeros are isolated if they are not present in the initial data. 

1 Introduction 

In a recent paper Sinai and Li [8] consider the initial-value problem for the 
three-dimensional incompressible Navier-Stokes equation, allowing the ve- 
locity field and the pressure to be complex- valued. They prove that, in this 
setting, there exist well-behaved (complex-valued) initial data for which the 
solution blows up in finite time. In this note we consider a similar problem 
for the viscous ID Burgers equation 

u t + uu x = u xx (1.1) 
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in R x (0, oo) with initial condition u(x, 0) = Uq(x), where we allow w to be 
complex-valued. A well-known fact about equation f 1 1.1 1) is that the transfor- 
mation u = —2v x /v, called the Cole-Hopf transformation, leads to standard 
heat equation v t = v xx for v. The singularities of u correspond to the zeros 
of v. For real valued functions, v cannot have zeros if they are not present in 
v(x, 0) and one sees immediately that for uq real and "sufficiently regular" the 
initial value problem for equation fll.ip has a unique smooth global solution 
(in some natural classes of functions), see [5]. This can, of course, also be 
seen without the use of the Cole-Hopf transformation, in a number of ways, 
since the equation fll.ip has a maximum principle and an energy estimate 
with respect to which it is subcritical. (The non-trivial scaling invariance of 
equation fll.ip is the same as for Navier-Stokes: u(x, t) — > \u(\x, X 2 t).) 

The maximum principle and the energy estimates are lost when we pass 
to complex-valued functions. At the same time, existence proofs based on 
perturbation theory and Picard iteration, such as in [6] or [7J, work also in 
the complex case, and it is therefore natural to expect that the proofs of 
local well-posedness for Navier-Stokes in critical (i. e. scale- invariant) spaces 
work also for equation fll.ip . without using its "complete integrability" . One 
can therefore expect local well-posedness of complex-valued equation f 1 1.1 1) 
in L 1 (by analogy with [6]) and, in fact, in (BMO)^ 1 (by analogy with [7]). 
With the Cole-Hopf transformation, the local L 1 well-posedness becomes 
completely transparent (see below). It it not quite so with the local (BMO)^ 1 
well-posedness, which shows the subtle nature of the well-posedness result in 
[7J. Our focus will be on global well-posedness, and therefore we will work 
with the L 1 space which is very simple and - as we will see - completely 
adequate for the problems we will consider. 

Since the zeros of v produce singularities of u, it is easy to find compactly 
supported smooth (complex-valued) initial data Uq for which the solution 
of equation fll.ip blows up in finite time, see Proposition 12.21 The conti- 
nuity argument used in the proof of this proposition allows us to formulate 
more general sufficient conditions on uq for the solution u to blow up, see 
Remark I2.3( i). On the other hand, in Proposition 12.11 we give a sufficient 
condition on u for the solution to converge to zero. Using these results, we 
can then explicitly describe the boundary, in some subsets of the space of 
initial data L 1 , between the basin of attraction of zero and the region from 
which the solutions blow-up, see Remark I2.5( iii). The behavior of solutions 
with the initial conditions on this boundary is then naturally of interest. 
These solutions are global and bounded and we describe their asymptotics 
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as t — > oo, see Proposition 12.41 and Remark I2.5( iii). 

Once we know that a solution can develop singularities, we can ask about 
the nature of the singular set. We will prove that, roughly speaking, if there 
are no singularities present in the initial data, then the set of singularities 
of the function u defined by the Cole-Hopf transformation from v (that is, 
u = —2v x /v) is always discrete in R x (0, oo). This follows from a theorem 
about zeros of complex- valued solutions of Id heat equation (Theorem 13. 3p . 
In a "typical situation" the number of singularities of such a solution u will 
be finite. However, as we show in Section [5j certain regular initial data yield 
solutions with infinitely many (isolated) singularities. We will also briefly 
address the question what "right-hand side" the singularities produce in a 
suitable weak formulation of the equation (see Section 0J. 

The solution of equation (11 .11) defined by u = —2v x /v is analytic outside a 
discrete set. This is not a typical behavior of solutions of non-linear parabolic 
equations with singularities. In fact, it is reasonable to expect that, for many 
equations, analyticity in the time variable will be destroyed in the whole time 
level {(x,t),t = t } if we have a singularity at time t at some point x . This 
conjecture is based on the study of singularities of the Complex Ginzburg- 
Landau equation in [9]. As far as we know, the issue has not been much 
studied. 

In the case of the dispersive regularization of Burgers equation, which 
is the KdV equation u t + uu x = u xxx , the singularities for complex- valued 
solutions are studied in [2J . Viscous Burgers equation with complex viscosity 
is studied by means of the Cole-Hopf transformation in [10] and [11] for 
a particular real initial condition, with the main focus on the behavior of 
singularities arising in complex time. 



2 Cole-Hopf transformation and singularities 

For a complex- valued u G L X (R) we define U(x) = f_ u(£,)d£ and v = 
exp(— U/2). Vice- versa, given a complex- valued v G W ' (R) (the space of 
all absolutely continuous functions that have the derivative in L 1 (1R)) with 
v(x) 7^ in R and v(x) — » 1 as x — > — oo, we let u = —2v x /v. For time- 
dependent functions on R we apply the above transformations at each time 
level. 

A well known simple calculation shows that u satisfies equation (II. ip if 
and only if v satisfies the standard heat equation vt = v xx , see for example 
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[5]. (If one does not impose the normalization v(x) -> 1 as x -> — oo, the 
function v is only determined up to a multiplicative factor depending on 
time, and the heat equation for v needs an extra term which would account 
for this, see [5].) 

We can now solve the initial value problem for equation (11.11) with a 
complex valued u G L 1 (M) as follows. Set v (x) = exp{— \J_ v>o(£)d£} 
and let v be the bounded solution of the heat equation with initial data 
vq. It is easy to check that there is T > such that \v(x, t)\ > e > in 
R x (0, T) and hence u = —2v x /v is a well-defined local-in-time solution of 
equation (II. ip with u(x, 0) = uq(x). 

Proposition 2.1. In the notation above, assume that Uq G L 1 (R) with 
J R |Im-uo| < 2tt. Then equation (II. ip has a global smooth solution u with 
u(x, 0) = u (x). If in addition \ f R lmuo\ < 2n, then swp x \u(x, t)\ — > as 
t — > oo. 

Proof. When J R |ImM | < 2ti, the function vq(x) = exp{— ~ J_ u (^)d^} 
takes values in a convex sector of the form {z, a < arg(^) < (3} with (3— a < n 
and it has finite nonzero limits as x — > ±oo. Thus for suitable real 9 we 
have Re (e i9 Vo(x)) > 0, x G R. Applying the strong maximum principle 
to Re (e ld v) (which solves the heat equation), we see that Re (e t9 v) > 0. 
Therefore v cannot vanish at any point (x,t) G R x (0, oo), proving the first 
statement. If | L Imit | < 2n, we can choose 9 so that Re(e* e f (±oo)) > E\ > 
0, and therefore for all large t we will have \v(x, t)\ > S\/2. Moreover, since 
v x also solves the heat equation and v 0>x G iv 1 (R), we have sup x \v x (x, t)\ — > 
as if: — > oo. These properties imply the second statement. □ 

Proposition 2.2. For each 5 > there exists a smooth, compactly sup- 
ported (complex-valued) uq with J R \uo\ < 2ir + 5 such that the solution of 
equation (11.11) with initial condition Uq blows up in finite time. 

Proof. We choose a smooth compactly supported non-negative ip with J R <p = 
2n + 5/2. Set u = —iip and let v be the Cole-Hopf transformation of u . 
The function v o satisfies: 

• v (x) = 1 for large negative x, 

• v (x) = exp(i(7r + 5/4)) for large positive x, 

• < a,rg(v (x)) < 7r + 5/4 for x G R. 
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Let v be the solution of the heat equation with initial data Vq. For 6 G R the 
function e~ ld v also solves the heat equation and choosing 9 > small enough 
we achieve that 

lim lm(e~ ie v (x)) < 0. 

x— >±oo 

It then follows that for a sufficiently large to > we have lm(e~ t0 v (x, to)) < 
for all i£8, Since the limits of e~ td v(x, t) as x — > ±oo are independent of t, 
comparing the trajectories of x i— ► e~ t9 v (x, t) for £ = and t = t we conclude 
that i> has to vanish at some point (xi,ti) with t\ G (0, to)- Consequently w 
has a singularity at (xi,ti). □ 

Remarks 2.3. (i) The above continuity argument can also be used to show 
that u develops a singularity whenever w G L 1 (R) satisfies | Llm«o| > 2tt 
and J R ImMo is not of the form 2it + 4/c7r, where k is an integer. 

(ii) It is perhaps worth pointing out a very non-local behavior implied by 
Propositions I2. II l2T2l Consider a compactly supported uq with Im« > and 
j R Im Mo = 7r + 5 for some small 5 > 0. Then the solution of equation (II. II) 
with initial condition uq exists for all time and converges to zero. Consider 
now the initial condition Uq(x) = u (x — a) +u (x + a). With initial condition 
Uq, the solution of (II. ip will blow up, no matter how large a is. If we take a 
very large, the solution will become very small in L°° before it starts growing 
again and blows up. (In fact, it is not hard to see that one can replace L°° 
by L p for a fixed p > 1 in the last sentence.) 

Proposition 2.4. Assumeuo G L 1 (R) is compactly supported, with | J M Im«o| 
j R | Imwol = 2tt. Then there exist a real y a and a complex (3 with \mf3 ^ 
such that the solution u of equation (II. ip with u(x, 0) = u (x) satisfies 

u(x, t) = ~ 2 + 0{\ ), (t^oc) (2.1) 

(x - y a V2t) + (3 \Jt 

uniformly in x G R. 

Proof. By Proposition ^. 11 the solution u is global. Let U (x) = u (£) d£ 
and v (x) = exp(— Uq(x)/2). Let L > be such that u vanishes outside 
[—L,L]. We have v (x) = 1 on (— oo, — L) and v (x) = — exp(— 1/2) on 
(L,oo), where I = f R Reu . We set F(x) = {2-k)- 1 ' 2 f* exp(-£ 2 /2) d£. 
(In terms of the erf function used in probability we can write F(x) = 
erf (x) — 1/2.) We note that the function F(^=) solves the heat equation in 

R x (0, oo) with the initial data \ sign (x). The fundamental solution T(x, t) 
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of the heat equation can be written as T(x, t) = We write the 

initial condition vq in the form 

v (x) = -a s 3^1 + b + Wo ( x ) (2.2) 

where a, b are chosen so that wq be compactly supported. This gives a = 
l+exp(— J/2) and b = (1— exp(— J/2))/2. Also note that Imu>o is continuous 
and does not change sign. Let t) be the solution of the heat equation in 
K. x (0, 00) with initial condition wq. The solution with the initial condition 
Vo is then 

v(x,t) = -aF(-^=)+b + w(x,t). (2.3) 



From the representation formula w(x, t) = Lw (y)T(x — y, t) dy we see that 
\w{x,t)\ + \Tt\w x {x,t)\ = 0{-=), (t^oo) (2.4) 

vt 

uniformly in x. To get further estimates for w in a simple way, we will use 
Appell's transformation and write 

w{x,t)=T(x,t)w{^,--), (2.5) 



where w is a function of x — x/t and t = —1/t (defined for (x,t) G 
(—00, 0)) which again satisfies the heat equation in x and t. We have 



x 



/ ~\ f \F(x — y,t) f xy ty 2 

w{x,t) = / w (y) dy = / w {y)exp{ — — + — ) dy , (2.6) 

Jr 1 1-''- t ) Jr - 1 

which shows that w can be analytically extended to MxR. Letting c = f R w , 
we see from (12.61) that w(0, 0) = c. Setting z = w — c, we can write 

v{x, t) = -a F{^=) + b + cr(x, t) + z{- -\) F{x, t) , (2.7) 
y It tt 



where z is smooth in R x (—00, 0] and z(0, 0) = 0. We now set a = b/a and 
(3 = - c /a. Observe that a G (-1/2, 1/2) and Imp ^ 0. From (ETj) we get 



n = -2^ = -2^^ ^ 2_L, (2.8) 



where F is evaluated at x/yM. Given any fixed R > and t < we can 
see that for \x\ < R, i > t we have \z$\ < ci, |^| < C2|x| + C3 1 1 1 . Together 
with (|2.8p . and after taking into account that Im/3 7^ 0, this gives 

1 + 0(4) -2 1 
n = _2^ 1^1 ^ = + 0(4=), (t->oo) (2.9) 

uniformly in regions {(x, t); < R}, where F is again evaluated at 

x/V2t. 

Let y a be the unique root of the equation F(y) = a and let us fix a (small) 
5i > 0. We note that when \x/y/2t—y a \ > 8\, then \aF{x/y/2b) — b\ > E\ > 0, 



and from (I2.3p . H2.4[) we see that, in the region {(x,t); \x/y2t — y a \ > Si}, 
one has 

« = -2- = 0(i=), (t-oo), 

uniformly in x (for (x,t) in the above region). Taking into account (12.91) . we 
see that it only remains to show that 

+ 0(^=), (*-»oo) (2.10) 



uniformly in {(2, t); \x/y2t—y a \ <5i}. (As above, F is evaluated at x/y2t.) 
For y e R we set 

F(y)-F{y a ) 
[V) F'(y)(y-y a ) 
with the understanding that K(y a ) = 1. Since T(x, t) = we nave 

—2 —2 —2 

3g^ = ^|^^ (2 ' U) 

Clearly 

K(2/) = l + 0(|y-yo|), (y^y a ). (2.12) 

Using the elementary formula 

^ ds 



k£ + P e + Jo (m-s) + sK}+(3y 
we see that (l2TT0ft follows from fl23TD and fl2~T2l . □ 
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Remarks 2.5. (i) We note that, with the notation used in the proof, we have 
Im/3 = - J R lmw . Under the assumptions of Proposition 12.41 the function 
Imtoo does not change sign and is integrable. If we do not assume that uq 
is compactly supported, it can easily happen that | J R lmwo\ = +00. We 
then have | Im/3| = +00 and in view of (12. ip it is natural to expect that in 
that case u(x,t) — > as t — > 00 uniformly in x. This can indeed be proved. 
If Uq is not compactly supported, but Jlmw is finite, we expect that the 
asymptotics of u(x,t) will be similar to (12.11) . with perhaps a slower rate of 
convergence. 

(ii) The constant y a is given by the equation F(y a ) = |tanh(|), with 
/ = f R Reu . In particular, y a = if and only if J R Rew = 0. We note 
that ^=4 is a steady-state solution of equation ( 11.11) . For / ^ we have 

^(y a ^/2t) = ^= 7^ 0, and we can interpret formula ( 12. ip as an "almost 
steady state solution", which is slowly drifting to ±00, at speed 

(iii) Consider a complex valued L 1 function uq supported in [-L, L], with 
J R I Initiol = / R I m ^o = 27T. Let G(u ,L,e) = {u E L 1 (1R) : \\u — uo\\ L i < 
e, uq is supported in [— L, L]}. From the above one can see that for suffi- 
ciently small e, one has, in the set 0(Uo, L, e), an explicit description of 
the boundary between the basin of attraction of the zero solution of equa- 
tion (11.11) and the region from which the solutions of equation ( II .ip blow 
up is finite time: The boundary (in O(u ,L,e)) is given by the equation 
/ R ImMo = 2n. (To be precise, for the proof of this one needs to augment 
the above propositions by a slightly modified version of Proposition 12.11 i n 
which we assume that uq is in O(uo, L,e), J R lmuQ < 2tt, and we allow 
L|Imw | < 2n + 5, where 5 = 5(L,e,-u ) > is sufficiently small. The 
restriction on the support of u Q is crucial in this step. We leave the details 
to the reader.) It is not hard to check that the large-time asymptotics of the 
solutions starting at the boundary (in O(u , L, e)) of the basin of attraction 
of the zero solution is given by the solutions described in Proposition 12.41 If 
we replace 0(uo,L,e) by (9(-u , 00, e) (i. e. we remove the restriction on the 
support of the perturbed function), the situation changes and the boundary 
is no longer described in a simple way. In addition, even when e is small, we 
expect that for O(uo, 00, e) some solutions at the boundary of the basin of at- 
traction of zero will have more complicated behavior, such as slow oscillations 
with large amplitude. 
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3 Nodal sets of caloric functions 

Let u be a bounded real-valued nontrivial solution of the heat equation u t = 
u xx in R x (0, oo). We define 

Z = {(x,t) G R x (0,oo); u(x,t) = 0} , 

Z rcg = {(x, t) G Z; u 2 + u 2 x ^0} , and 

•^sing ^ \ ^rcg ■ 

The analyticity of u implies that Z sing is discrete and that Z is locally a 
regular real analytic curve in a neighborhood of each point (x ,t ) in Z reg . 
By a regular (real) analytic curve C in an open set U Clx (0, oo) we mean 
a one-dimensional analytic (imbedded) submanifold of U with C\C C cM7. 

Lemma 3.1. In t/ie notation introduced above, the regular analytic curves 
describing Z in a neighborhood of(xo,to) G Z rcg can be analytically continued 
through the points of Z sing . In other words, Z is a (locally finite) union of 
regular analytic curves m R x (0, oo). 

Proof. We first recall some facts about caloric polynomials. As usual in the 
parabolic setting we say that a function f(x, t) is parabolically m— homogeneous 
if f(\x,\ 2 t) = \ m f(x,t) for A > 0. The m— th caloric polynomial is a 
parabolically m— homogeneous polynomial satisfying the heat equation. It is 
unique, modulo a multiplicative factor, and can be given for example by 

[m/2] ( 

(3-D 

The polynomial P m (x, —1) is the m— th Hermite polynomial. For m even, 
m = 2k, the polynomial P m is of the form 

P m (x, t) = (x 2 + ait) • • • (x 2 + a k t), (3.2) 

with < a\ < ■ ■ ■ < afc. For m odd, m — 2k + 1, P„ t is of the form 

P m (x, t) = x(x 2 + ait) • • • (x 2 + a^), (3.3) 

with < ai < • ■ • < afe. (The a^-'s may be different for different m, of course.) 

In a neighborhood of a point (x , t ) G ^ s in g we can write u as a convergent 
series 

w(x,t) = a m P rn (x -x ,t- t ) + a m+ iP m+ i(x - x ,t - t ) H (3.4) 
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where m > 3, a m ^ 0. In what follows we will assume that m is odd, 
m = 2k + l. (The proof for m even is similar and, in fact, easier.) We change 
coordinates so that (xq, to) corresponds to (0, 0) in the new coordinates, which 
we still denote (x,t). We let t = —y 2 ,ai = b\, ■ ■ • , = b 2 ., bj > 0. The 
equation u(x,t) = can be written as 

u(x, -y 2 ) = x(x - hy)(x + b x y) ■ ■ ■ (x - b k y)(x + b k y) + R(x, y) = 0, (3.5) 

where R(x, y) is analytic in a neighborhood of (0, 0) with vanishing deriva- 
tives of order 1,2, ...m. Letting 6 — 0, we will look for analytic curves 
x = x(y) of the form x(y) = by + y 2 f(y) (with b = ±bj, j = 0,1, ... ,k) 
defined for small y on which u(x, —y 2 ) vanishes. Substituting the expression 
x(y) = by + y 2 f(y) in equation (I3.5p . it is easy to check that we get an 
equation of the form 

f(y) = F(y,f(y)) (3.6) 

where F = F(y, /) is analytic in / and depends on / only through yf. 
Therefore fo := F(0,f) is independent of / and §j(0, /) = 0. Applying 
the standard implicit function theorem one shows that equation (I3.6P has an 
analytic solution / defined on a neighborhood of with /(0) = Jo- 
Observe that with y ^ fixed, the curves found above give m different 
solutions of (13.51) . Therefore, by the Weierstrass preparation theorem, they 
yield all solutions of ( 13. 5ft in a neighborhood of (0, 0). (One can also use the 
Malgrange preparation theorem.) 

To finish the proof, we note that for j > 1, instead of writing x = x(y) 
we can write y = y(x) and the equation t = —y 2 = —(y(x)) 2 then defines the 
analytic branch of Z which has contact of the second order with the parabola 
t = —x 2 /cij. When j = we note that the function in equation ( 13.61) is of 
the form F(y, f) = F(—y 2 , f), and hence the corresponding curve is of the 
form x — x(t) — tf(t), with an analytic / satisfying /(0) 7^ 0. □ 

Remark 3.2. It is clear that the proof of the lemma also works when the 
equation has lower-order terms and analytic coefficients. Although we did 
not find the precise statement of the lemma in the literature, we assume 
it is known to experts. For example, it follows easily from the analysis of 
nodal sets in [T], where the method of Newton polygons is used. Also, once 
the specific form of the caloric polynomials is taken into account, the lemma 
can be derived easily from general principles used in algebraic geometry for 
"desingularization" . Nevertheless, we think that the elementary proof above 
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is still of some interest and we have included it for completeness. We remark 
that even if we allow nonanalytic variable coefficients, Z is still a finite union 
of regular C 1 curves in a neighborhood of any point in Z s { ng , see [3]. 

Theorem 3.3. Let v be a bounded complex-valued solution of the heat equa- 
tion mix (0, oo). Assume v has no zeros in some neighborhood o/R x {0}. 
Then all zeros ofvinWx (0, oo) are isolated. 

Proof. Let v = Vi+iv 2 , Z\ = {vi = 0}, Z 2 = {v 2 = 0}. Assume Z 1 C\Z 2 has 
an accumulation point (xo, to) inside Rx (0, oo). By Lemma [37TI we know that 
Z\ is a locally finite union of regular analytic curves in R x (0, oo). Consider 
the curves passing through (xo,to)- Clearly v 2 has infinitely many zeros 
accumulating at (xq, to) on one of the curves, let's call it C. By analyticity, 
v 2 vanishes on C. The curve C cannot be closed, for otherwise the maximum 
principle would imply that both V\ and v 2 vanish in the interior of C, which 
is impossible by our assumption and analyticity. Hence we can parametrize 
C by a parameter s G (—00,00). Also, since C is a regular analytic curve, 
we have C\Cclx {0}, thus our assumption implies that C\C = 0. Now 
either the time coordinate t has a strict local minimum on C or we can choose 
the parametrization so that t(s) is monotone nonincreasing for large s and 
x(s) approaches 00 or —00 as s — >• 00. In either case, we find a (bounded or 
unbounded) domain in R x (0, 00) such that both functions v% and v 2 vanish 
on its parabolic boundary. Since they are bounded, the maximum principle 
[1] implies that they vanish on a nonempty open set, hence on R x (0, 00), 
and we again have a contradiction to our assumption. 

□ 

It is clear that the proof of Theorem 13.31 works without much change also 
for complex-valued harmonic functions in a half-plane. 

4 Additional comments on the singularities 

Given complex- valued initial data uq G L x (R) for equation 11.11 and con- 
structing the solution u of the initial-value problem by means of the Cole- 
Hopf transformation as in Section 2 by setting u = —2v x /v, we see from 
Theorem 13.31 that the singularities of u are isolated. It is natural to ask if 
equation (11.11) is satisfied in some weak sense across the singularities, or if 
the singularities introduce a non-trivial "right-hand side", i. e. we want to 
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calculate the distribution / given by u t + uu x — u xx = f, where the left-hand 
side requires a suitable interpretation. Clearly / should be supported in the 
singular set. Even if we write the operator u t + uu x — u xx as u t + (u 2 /2) x — u xx 
the definition of / can still be somewhat ambiguous, since u and u 2 are not 
locally integrable in a neighborhood of a singularity. We suggest one possible 
interpretation. For simplicity we will consider only the simplest case when 
the function v defining u has a simple zero at the singularity (x ,to), i. e. 
v(x, t) = a(x — xo) + b(t — to) + 0((x — xq) 2 + (t — to) 2 ), with a, b complex and 
linearly independent over R. We consider a smooth test function ip = ip(x, t) 
supported in a small neighborhood of (xo,t ), so that no other singularity is 
present in the support of if. We want to define 

1=1 (-uip t - u 2 ip x /2 - uip xx ) dxdt. (4.1) 

For t 7^ to we let 

h(t) = / (—u(x,t)(p t (x,t) - u 2 (x,t)ip x (x,t)/2 - u(x,t)<p xx (x,t)) dx (4.2) 



Using well-known facts about the behavior of the distributions l/(x ± is) 
and l/(x ± is) 2 as e 0, together with a change of variables v(x, t) = y (for 
a fixed t) one can see that h(t) has one-sided limits as t — > to- Therefore it 
seems to be natural to define integral / in expression ( 14.11) as 

I = h(t) dt = lim ( / h(t) dt+ h{t) dt.) (4.3) 

J0 Jo Jto+T 

Integration by parts, the equation satisfied by u, and the specific form of the 
singularity of u now give 

I = lim / \u(x, t + r)(p(x, to + r) — u(x, to — r)ip(x, to — r)l dx = ±4Tri(p(xo), 
t ^°Jr 

(4.4) 

where the correct sign is the same as the sign of the imaginary part of ab. In 
other words, in a neighborhood of (xo,t ) we have, in some sense, 

u t + uu x - u xx = ±4iri8( XOjto ), (4.5) 

where we use the usual notation 5r XO)to \ for the Dirac distribution at the 
point (xo,to). Therefore we cannot interpret u as a global weak solution of 
equation (II. ip . 
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5 Infinitely many singularities 



If is not hard to show that, typically, the solution u will only have finitely 
many singularities. In fact, one can check easily that a sufficient condition for 
u to have only finitely many singularities is that J R Im«o(x) dx is not of the 
form 27T+4/c7r with k an integer. We now show that, on the other hand, there 
are solutions with regular initial data having infinitely many singularities. 
This is an immediate consequence of Proposition 15.11 below. We recall that 
we denote by W ' (R) the space of all functions on R which are absolutely 
continuous with the derivative in L 1 (M). (In particular, constant functions 
belong to W ' .) 

Proposition 5.1. There exists a smooth (complex-valued) function vq G 
Wq' (R) such that t> (— oo) = 1, \v (x)\ > £o > for any x G R, and the 
solution v of the heat equation with v(-,0) = Vq vanishes at infinitely many 
points (0,Tfc), with r k — > oo. 

Proof. First choose a smooth real- valued odd function w G W^'^R) such 
that w(— oo) = 1 and w > on (— oo, 0). The solution v\ of the heat equation 
with vi(-, 0) = w has a unique zero at x = for each t. We shall next find 
a smooth real-valued function z G W ' (R) such that z(±oo) = 0, z(0) ^ 0, 
and the solution t> 2 of the heat equation with f 2 (-, 0) = z vanishes at points 
(0, Tfc) with Tfc — > oo. From this the conclusion of the proposition follows 
upon setting v = w + iz. 

We first choose sequences Rk > 0, G (0, 1) with the following properties: 

( al ) J2T=l e k < OO, 

(a2) R k+1 > R k and R k+1 e k+1 > J2j=i e i( R i + l ) (A; = 1, 2 . . . ). 
Next, for each k we choose a smooth function z k such that 
(a3) z k = e k on [-R k , R k ], z k = on R \ [-R k - 1, R k + 1), 
(a4) < z k < e k and \z' k \ < 2e k on R. 

We will show that if x k is a suitably chosen sequence, then the function 

oo 

z{x) = Y,^) k+l ^{x-x k ) (5.1) 
k=l 
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has the desired properties. 

The sequence Xk will be constructed so that, in particular, 



Xj\ > + Rj + Rj-i + 2 



(5.2) 



for j = 2, 3 ... . This guarantees that the functions Zf-( ■ — Xk) have non- 
overlapping supports, hence, by (al), (a3), and (a4), z is a smooth function 
in Wo'^M) satisfying ^(±oo) = 0. In addition to (15. 2p . we need to ensure 
that the function 



has infinitely many sign changes. 

We shall recursively construct sequences x^, d k such that (— l) k+1 x k 
tk, dk are nonnegative and increasing with k and the following statement is 
satisfied for each k: 

(Ik) Relations (15. 2p hold for j — 2, . . . , k and, with any choice of Xk+i, Xk+2, ■ ■ ■ 
satisfying (15.21) for j = k + 1, k + 2, . . . and \xk+i\ > dk, one has 



Take x\ — 0, t\ — 1. It is obvious that (Ii) is satisfied if d\ is sufficiently 
large. We fix such a di satisfying also d\ > Ri + R 2 + 2. 

Assume that Xj, tj, dj have been constructed for j — 1, . . . , k. To define 
the next terms assume for definiteness that k is odd (in case it is even, the 
construction is analogous). Set Xk+i = —dk- Assuming Xk+2, 2^+3, • • • , are 
any numbers satisfying ( 15.21) for j = k + 2, k + 2, . . . and \xk+2\ > dk+i, with 
d k+ i to be specified below, we use (a3), (a4) to estimate 




(5.3) 






(5.4) 



(5.5) 
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In view of (a2), if t = t k+1 > t k is large enough, the expression in (15.41) 
is negative. Fixing such a t k+ i and subsequently choosing a large enough 
dk+i > dk + Rk+i + Rk + 2, we make the whole expression in (15.41) . (15.51) 
negative. Hence (Ifc+i) is satisfied. 

It is obvious that with sequences x k and t k resulting from the above 
construction, the function z has all the desired properties. In particular, 
v 2 (0,t) has a zero r k in (t k ,t k+ i) for each k. □ 

Remarks 5.2. (i) It is not difficult to check that with the initial data z 
constructed above, the solution t> 2 of the heat equation has a unique zero 
x(t) for each t > (and x(f) changes sign infinitely many times as t — > oo). 

(ii) It is conceivable that if the initial condition w is compactly supported, 
the solution u of equation (II. ip given by the Cole-Hopf transformation will 
always have only finitely many singularities. For example, the construction 
above cannot be carried out if we demand that z be compactly supported. 
This can be seen from the following observation: If z solves the heat equation 
in M x (0, oo) with compactly supported initial data and t — > z(0,t) has 
infinitely many zeros, then z(0,t) = for all t > 0. To see this we recall 
that Appell's transformation z of z is defined in R x (— oo,0) by z(x,t) = 
T(x, t)z(x/t, —l/t). One can see from formula (12.61) applied to z that z has an 
analytic extension tolxR and the assumptions on z imply that t — > z(0, t) 
has infinitely many roots accumulating at 0. Hence z(0, t) = for all t, which 
implies our statement. 

Acknowledgments. The authors would like to thank Luis Escauriaza for 
a valuable discussion. 
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